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Raman induced soliton self-frequency shift in microresonator Kerr frequency combs 
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The formation of temporal dissipative solitons in continuous wave laser driven microresonators 
enables the generation of coherent, broadband and spectrally smooth optical frequency combs as 
well as femtosecond pulses with compact form factor. Here we report for the first time on the 
observation of a Raman-induced soliton self-frequency shift for a microresonator soliton. The Raman 
effect manifests itself in amorphous SiN microresonator based single soliton states by a spectrum 
that is hyperbolic secant in shape, but whose center is spectrally red-shifted (i.e. offset) from the 
continuous wave pump laser. The shift is theoretically described by the first order shock term of the 
material’s Raman response, and we infer a Raman shock time of 20 fs for amorphous SiN. Moreover, 
we observe that the Raman induced frequency shift can lead to a cancellation or overcompensation 
of the soliton recoil caused by the formation of a (coherent) dispersive wave. The observations are in 
excellent agreement with numerical simulations based on the Lugiato-Lefever equation (LLE) with 
a Raman shock term. Our results contribute to the understanding of Kerr frequency combs in the 
soliton regime, enable to substantially improve the accuracy of modeling and are relevant to the 
fundamental timing jitter of microresonator solitons. 


Introduction. — Microresonator based optical fre¬ 
quency combs (Kerr combs) [1, 2] enable optical fre¬ 
quency comb generation from a continuous wave (GW) 
laser, with repetition rates in the microwave domain (>10 
GHz), and broad spectral bandwidth [3-5]. Recently, a 
qualitatively new operation regime has been discovered 
[6], in which the parametrically generated comb seeds 
the formation of a temporal dissipative (cavity) soliton 
[7-9]. Such temporal dissipative solitons, first externally 
induced in fiber cavities [10], have been observed in Kerr 
frequency comb experiments using crystalline microres¬ 
onators [6, 11] and have recently also been generated in 
photonic chip based silicon nitride (SiN) integrated res¬ 
onators [5, 12] . Soliton based microresonator frequency 
combs have several attractive features, in particular be¬ 
ing fully coherent, having smooth envelopes and giving 
access to femtosecond pulses. Indeed, the short dura¬ 
tion of temporal solitons in crystalline microresonators 
have been used for external fiber broadening [6, 13] and 
have allowed to achieve self-referencing (i.e. determina¬ 
tion of the comb’s carrier envelope frequency) enabling 
to count the cycles of light [13[. Moreover, taking ad¬ 
vantage of dispersion engineering [3, 14], the presence of 
third (and higher) order dispersion allows (coherent) dis¬ 
persive waves (DWs) to be generated [12] via the effect 
of soliton induced Gherenkov radiation [15] . This process 
has been used to create a photonic chip based (coherent) 
frequency comb in a SiN microresonator spanning 2/3 
of an octave at electronically detectable mode spacing. 
However, while advances in theoretical simulations of the 
soliton regime have occured [14, 16, 17] a key underlying 
question is what other effects play a role in soliton for¬ 
mation and need to be included into the simulations. 

In contrast to “Turing rolls” [18, 19] dominating combs 
in the modulation instability (MI) regime, cavity soli¬ 
tons carry intense peak power and ultrashort duration 


such that in principle it could excite higher-order non¬ 
linear effects such as the self-steepening effect and Ra¬ 
man scattering. The latter usually exists and is broad¬ 
band in nature in amorphous oxide and nitride mate¬ 
rials such as silica and SiN, which relates to the vi¬ 
brational material response of the cubic nonlinearity. 
Material Raman studies revealed typical active modes 
of e.g. mono-silicon bonds (520 cm“^), silicon-oxygen 
bonds (450 cm“^, broadband) [20] and silicon-nitrogen 
bonds (400 cm“^ in amorphous Si 3 N 4 ) [21] . In silica 
fibers, the interplay of the Raman effect with an ul¬ 
trashort femtosecond pulse has been widely investigated 
[22, 23] in the context of fiber-based supercontinuum gen¬ 
eration and soliton propagation [24] . The observed phe¬ 
nomena in soliton propagation in fibers include among 
others soliton self-frequency shift that continuously shifts 
the whole pulse spectrum to the red upon propagation 
[25, 26] , soliton fission when a high-order soliton is split 
into fundamental solitons [27] , and Raman induced pulse 
compressions in dispersive media [28] . While well studied 
in fiber, to date however the role of the Raman effect in 
microresonator soliton formation has not been observed. 
Soliton formation in crystalline materials has not shown 
any evidence [6, 11[. Recently, for amorphous silica- and 
SiN-based microresonators, the Raman effect was esti¬ 
mated by numerical simulations [29-31], and it was found 
that a comb spectral shift can occur when a Raman- 
nonlinearity assisted single temporal dissipative soliton 
is formed. In this letter, we report for the first time to 
the authors’ best knowledge the direct experimental ob¬ 
servation of Raman induced soliton self-frequency shift 
in microresonator solitons. 

Experiments. — The study is based on SiN mi¬ 
croresonators [32] in which recently temporal solitons 
were generated [12] . The first set of samples was fab¬ 
ricated by a newly developed “photonic damascene pro- 
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FIG. 1. (a) Diagram of the experimental setup. AFG, arbi¬ 
trary function generator; ECDL, external cavity diode laser; 
EDFA, erbium doped fiber amplifier; EPC, fiber polariza¬ 
tion controller; OSA, optical spectrum analyzer; OSG, os¬ 
cilloscope; PD, photodiode; (b) Optical microscope image of 
the SiN microresonator, with a ca. 100 GHz mode spacing, 
(c) Illustration of the pump laser detuning excitation scheme 
(reducing the pump frequency). The gray line shows the Kerr 
nonlinearity tilted cavity resonance prohle with bistability re¬ 
ferring to the GW, the blue line shows the trace of intracavity 
peak power with an increase of the cavity resonance-pump de¬ 
tuning S {S = fo — fp , where fo and /p are the resonance and 
pump frequency), the red line shows the single soliton exis¬ 
tence range, (d) Illustration of the Raman induced spectral 
red-shift of the soliton frequency comb (in red line) compared 
to the case without Raman contributions (blue dashed). 


cess” [33]. Samples have a free spectral range (FSR) 
of 75-100 GHz on the fundamental TEqo mode fam¬ 
ily. Resonators have nominal SiN core height of 0.9 ^m 
and width 1.65 /im. The resonance linewidth ^ around 
1550 nm is ~ 200 MHz that corresponds to a loaded 
Q-factor of ~ 10®. The dispersion landscape is given 
by - (wo + Di ■ n) = i e N 

[34], where /i indicates relative mode index to the 
pumped mode = 0), FSR = ^, is the resonance 
frequency of the resonator, and |;i=o- At around 

1550 nm, parameters = 1 — 2 MHz (anomalous group 
velocity dispersion) and = ^(1 kHz) are character¬ 
ized. 

Fig. 1(a) shows the experiment setup. The CW pump 
is red-tuned to sweep over a resonance (at 1550 nm) 
with constant speed using an AFG (1 nm/s). The pump 
power is 1-3 W on chip. Usually, transitions to soliton 
states can be identified from the generated comb light 
[6] that shows step-like pattern, indicating the forma¬ 
tion of temporal dissipative solitons (e.g. Fig. 1(c)). 
Multiple solitons and single soliton are deterministically 
and repeatedly achievable by properly stopping the AFG- 


tuning. Typical single soliton based frequency combs 
(soliton combs) in our SiN microresonators are shown 
in Fig. 2(a-c). The comb spectral envelope shows sym¬ 
metric and smooth sech^-shape with a 3-dB bandwidth 
of ^ 5.5 THz, despite minor irregularities (e.g. the spike 
at 198.5 THz) caused by mode crossings [11[. One strik¬ 
ing difference to prior work in crystalline resonators [6] is 
a global red-shift of the solitonic part of the comb spec¬ 
trum, compared to the GW pump frequency, as shown in 
Fig. 2(a). 

Once the single dissipative soliton is formed, one can 
further explore it’s laser detuning dependence, by tuning 
the pump either to red or blue wavelengths, while the dis¬ 
sipative soliton still persists. The range of laser detuning 
over which the soliton state exists (i.e. the “soliton step 
length”) usually corresponds to several GHz. 

We then set the blue end of the single soliton step 
as the initial detuning (5^) and tune the pump wave¬ 
length over the entire laser detuning step length. With 
an increase of the detuning, the slowly evolving comb 
shows two apparent trends: first, the single soliton comb 
spectrum is broadened as the soliton duration is short¬ 
ened and second, the global comb spectral red-shift is 
increased. Figure 2(a, b) shows two comb spectra (in lin¬ 
ear scale) at two positions within the single soliton step. 
The comb red-shift in (a) is larger than that in (b) since 
the laser detuning S is increased. Figure 2(c) shows the 
overall spectral broadening (in logarithmic scale) upon 
an increase of S. For comparison, a high noise MI comb 
is also measured before the pump wavelength is tuned 
into the soliton existence range, see Fig. 2(d), in which 
the spectral shift is not observed (in agreement with the 
understanding that this regime does not produce short 
pulses inside the cavity). Moreover, the spectral red-shift 
on single soliton combs is found to be changed by laser 
tuning within the single soliton step and such tunability 
is reversible by reversing the scan direction in contrast 
to the continuous red-shift in fiber optics [22, 23] . The 
comb red-shift as a function of the detuning is shown in 
Fig. 2(e) for different pump powers. The shift ranges 
from 0.5 THz (~ 5 x FSR) to 1.75 THz (~ 17.5 x FSR) 
and exhibits almost linear dependence with the laser de¬ 
tuning in the single soliton existence range. The soliton 
step length is also increased with an increase of the pump 
power, shown in Fig. 2(e). 

We attribute the comb spectral red-shift to the Raman 
induced soliton self-frequency shift [22, 23, 35], which 
in the studied SiN microresonator is well isolated and 
doesn’t mix with higher-order dispersion effects, since the 
resonator produces weak third order dispersion and effec¬ 
tively suppresses the DW generation via soliton induced 
Gherenkov radiation [15]. In addition, highly symmetric 
comb spectra imply weak self-steepening effects. 

As a third order nonlinear effect, the Raman induced 
soliton self-frequency shift is in principle scaled by the 
soliton peak power, which matches the trend we observed 
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FIG. 2. (a, b) Experimentally generated Kerr combs in the 
single temporal soliton regime in a 100 GHz mode spacing 
SiN microresonator, shown on a linear scale. The envelopes 
are from numerical simulations. The comb power (y-axis) is 
normalized to the GW pump power. For soliton comb sim¬ 
ulations, we used = 2 MHz, ^ = 4 kHz, ^ = 350 MHz, 
(5i = 2.1 GHz. (c) Measured soliton combs shown in logarith¬ 
mic scale, red comb (same to (a)) has <5 = di -|- 1 GHz, the 
green one corresponds to 5 = di -I- 0.5 GHz, and the blue one 
(same to (b)) corresponds S = 5i. (d) Experimental measure¬ 
ment and simulation (the envelope) of a high noise MI comb, 
with the detuning <5 = 0.35 GHz outside the soliton existence 
range, (e) Variation of the comb spectral red-shift as a func¬ 
tion of the detuning (ai-axis: 5 — 5i), over the single soliton 
step at three pump powers: 1.9 W, 2.3 W and 2.7 W. Note 
that the trend from the simulation is discrete (gray line) as 
we mark the specific comb line of highest power (despite the 
strong pump line). The spacing corresponds to the FSR. ^ 


in experiments. Analytically, the single temporal dissipa¬ 
tive soliton consists of a CW background solution T'cw 
and a superimposed temporal soliton 'Fg =‘^sech^(^) 

where the soliton peak power is ^ oc (n, n 2 , I4ff 

are waveguide refractive index (RI), material nonlinear 
RI and effective mode volume, respectively), i.e. linearly 


proportional to the laser detuning <5 [6, 8]. The soliton 
pulse duration, on the other hand, is At oc which 
is inverse proportional to the detuning [6]. Therefore, 
benefitting from the increase of the soliton power and 
the decrease of the pulse duration upon a red-side pump 
wavelength detuning, the enhancement of the comb red- 
shift and the spectral broadening are expected. 

Theory and numerical simulations — We numerically 
simulate the Kerr comb generation in SiN microres¬ 
onators, based on LLE [b, 7, 36, 37], but include the 
Raman effect [29-31]: 


dA{n,t) 


dt 




^ex'^inl^—0 ( T ^27r:5 ) A 


■ igF [(1 - /r)|A| 2A + /R(hR((^) 0 jA^A] (1) 


where |A(/r)p indicates the number of photons in the 
mode with index /r and t is the propagation axis (i.e. 
the slow time frame). marks the envelope written 
in frequency-like n domain. A{(j)) indicates the envelope 
in the phase 4> = domain, where r is the fast time 
frame and is the round-trip time. Moreover, the 

input drive is Sin = (Pi„/Ikuo)^^^ where Pin is the pump 
power. The cavity decay rate k = Kgx + kq is the sum of 
intrinsic decay rate kq and coupling rate to the waveguide 
Kex- The nonlinear Kerr coupling is g(fi) = 

Here, //-dependent parameters (n, n 2 and Ves) imply that 
the self-steepening effects can automatically be included 
in the model. F indicates the Fourier transform, /r is 
the Raman fraction and = Pi“^/iR(r) is the Ra¬ 

man response scaled in (j) domain where h^{T) is in the 
physical (fast) time (r) domain. Therefore, the Raman 
spectrum has /ir(//) = hR(P[~^a;), where uj marks the 
physical frequency domain. In the case that Raman ac¬ 
tive modes are sufficiently high in frequency and exceed 
the bandwidth of the soliton pulse, the resonant Raman 
effect plays a negligible role [38] . Then the Raman term 
can be simplified only containing the instantaneous re¬ 
sponse (direct current (DC) component) that contributes 
to the electronic Kerr effects, and the first-order compo¬ 
nent identified as the Raman shock term [38]: 

= ( 2 ) 

fR d(j) 

where ())r = = ih'j^ifJ, = 0) and tr is the material 

Raman shock time. For the studied SiN microresonator 
samples, the fundamental single soliton is formed with a 
3-dB bandwidth ^5.5 THz that is lower than the fre¬ 
quency of material Raman modes (~ 12 THz). 

Theoretically, the soliton self-frequency shift is fully 
described by the Raman shock term [25, 35[ that di¬ 
rectly affects the pulse phase, thereby inducing a tempo¬ 
ral phase shift that correspondingly red-shifts the pulse 
spectrum. It is essential to delineate this process from 
the self-steepening effect that also introduces a shock 
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FIG. 3. (a) Illustration of the cancellation effect between 

the Raman induced spectral red-shift and soliton blue-recoil 
through DW emissions on the long wavelength side. (b,c) Ex¬ 
perimental generation and simulation of single soliton combs 
with coherent DWs in two different SiN microresonator ge¬ 
ometries (both having a nominal SiN thickness of 0.8 pm 
and mode spacing of ~190 GHz). Panel (b) shows the sin¬ 
gle soliton comb from the sample with a SiN core width of 
1.8 pm (sample taken from [12]), the comb 3-dB bandwidth 
is ~ 10 THz, Dint for the simulation contains ^ = 2.2 MHz, 
^ = 18 kHz, ^ = —350 Hz, and the pumped resonance is 
at 1560 nm employing 1 W of power Panel (c) shows a soli¬ 
ton comb in a different sample with SiN core width of 2.0 pm 
and the 3-dB bandwith is ~ 8 THz. Simulated dispersion pa¬ 
rameters ^ = 3.2 MHz, ^ = 26 kHz, = -340 Hz. The 
pumped resonance is at 1554 nm and the power is 1 W. 
Dint and Vea are calculated by finite element modelling us¬ 
ing C0 MSDL(r), while D 2 is also measured [34], showing close 
agreement to simulations. Blue dashed lines indicate comb 
envelopes without the Raman contribution. 


term, but works directly on the pulse temporal ampli¬ 
tude. An explicit self-steepening term could be recog¬ 
nized when the nonlinear coupling is written in (j) domain, 
i.e. F-^[g] = s ■ go, where go = g{g, = 0). § = 1 - 
indicates the self-steepening operator, where the shock 
time (scaled in </) domain) is ~ 

dicates the first-order variation of I4g(p) (in p) around 

p = 0. 

We expand Eq.(l) into amplitude and phase dynamics 
of A{(j)), using the ansatz As¬ 

suming the self-steepening as well as the Raman effect 
only make perturbations to the cavity single soliton that 
is the eigen-solution of the standard LLE (Di >2 = 0, 


(/)s = 0, and /u = 0), we have following two relations: 
dAu , 2 ^^ r, r j. du 


where Au and A(f> indicate the amplitude and phase per¬ 
turbations. It is clear that while the phase perturbation 
is induced by the Raman shock term (()r, the amplitude 
perturbation is related to c/ig. Self steepening thereby 
gives rise to shock front on the soliton pulse that distorts 
the spectrum and makes it asymmetric (opposite to the 
Raman case). In experiments, a highly symmetric comb 
envelope (Fig. 2) actually implies weak self-steepening ef¬ 
fects. Reasons include that, first the comb bandwidth is 
small compared to the pump frequency, and second, the 
non-zero AV helps reducing the overall self-steepening 
effects. 

Figure 2 contains a set of numerical simulations which 
shows almost perfect agreement to the experiments. We 
set the Raman fraction /r = 20 %. Importantly, the Ra¬ 
man shock time tr in SiN is so far unknown and has not 
been reported to the authors’ best knowledge. However, 
close agremeent between theory and experiment and pre¬ 
cise knowledge of experimental parameters allow us to 
determine the parameter to be 27rrR ^ 20 fs. Moreover, 
as independent verifications, this value is found to agree 
between different SiN microresonators, including those 
with significant different structural geometries. By lin¬ 
early tuning S, the comb evolution is simulated and tran¬ 
sitions from primary comb, noisy MI comb to temporal 
dissipative soliton states are identified. The trace of the 
intracavity peak power (cf. Fig. 1(c)) shows the sin¬ 
gle soliton existence range 2TrS = 6 — 11k. Figure 2(a,b) 
show simulated comb envelopes in the soliton state, both 
with significant red-shifts. In a high noise MI comb [6] 
(Fig. 2(d)), however, the chaotic cavity light pattern is 
identified which is not as stable and intense as a cav¬ 
ity soliton, and the comb spectral red-shift is absent in 
agreement with experiment. Simulations also verified the 
trend of the comb red-shift as a function of the detuning 
S, see Fig. 2(e). 

It is interesting to note that the Raman shock time 
we extracted for SiN in this work (27rTR ~ 20 fs) is ac¬ 
tually much smaller than that in silica, the latter be¬ 
ing 27rrR « 27r X 89 fs [20]. A large shock time enables 
the observation of soliton self-frequency shift in silica- 
based fibers with short lengths (^(10) cm) and intense 
pump powers (^(1) kW) [39], while in SiN-based waveg¬ 
uides (not resonators), the Raman effect is not observed 
[40, 41]. However, in microresonators, due to the cavity 
building (finesse ^ « ^(100- 1000)), the light- 

material interaction length is effectively increased (to 
^(1) m) and meanwhile the pulse peak power is dra¬ 
matically promoted (e.g. in Fig. 1(c) the intracavity 
soliton peak power is up to ca. 3 kW under a CW pump 
of 2 W). Therefore, the weak Raman effect can be ex¬ 
cited. Moreover, the single soliton formation in a mi- 
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croresonator promises a persistent fundamental soliton 
pulse that could naturally suppress high-order Raman ef¬ 
fects (e.g. the soliton fission) and leave the pure effect of 
soliton self-frequency shift, despite that the comb band¬ 
width may reach the frequency of Raman active modes. 
To this point, the Raman effect actually assists the for¬ 
mation of the soliton comb, provided that the Raman 
gain doesn’t exceed the cavity loss (see a recent careful 
study in [29]). Therefore, even though the shock time is 
much shorter than that in silica glass, we could still ob¬ 
serve the Raman induced soliton self-frequency shift in 
SiN microresonators. 

Cancellation of the soliton spectral recoil by the Raman 
shift. — We next investigate how the Raman effect in¬ 
fluences cavity soliton with the presence of higher-order 
dispersion effects. For this we use two other sets of SiN 
microresonators with two different structural geometries 
that allow shorter soliton pulses and enable coherent DW 
generation, as recently demonstrated in [12]. Third order 
dispersion ^ = ^(10) kHz is characterized in the waveg¬ 
uide, on the fundamental TMqo mode. A zero dispersion 
wavelength is estimated at ca. 1700 nm. Experimentally 
pumped at 1550 nm a single soliton comb is generated 
in both samples, and more importantly, a coherent DW 
is identified at around 195 THz (cf. [12] for more exper¬ 
imental details). Since a coherent DW is spectrally sit¬ 
uated in the normal dispersion regime, the overall comb 
bandwidth is significantly extended to ~ 75 THz. The 
comparison of the two spectra in Fig. 3(b) and Fig. 3(c) 
shows the influence of dispersion on the position of the 
DW frequency. 

It is well understood that the emission of a DW will 
lead to the soliton spectral recoil [15, 42], such that the 
comb profile around the CW pump frequency will be 
shifted away from the DW, as illustrated schematically 
in Fig. 3(a). However, in above mentioned SiN microres¬ 
onators, the blue-recoil (resulting from the DW being on 
the red side) was surprisingly not observed [12], see Fig. 
3(b), and even a red-shift (Fig. 3(c)) of the comb pro¬ 
file is seen. The explanation for this unusual observation 
is that the Raman induced soliton self-frequency shift 
cancels and even overcompensates the soliton blue-recoil. 
Indeed, a similar effect has been observed in the case of 
soliton propagation in optical fibers [43] . We compare 
the two spectra to numerical simulations including both 
the Raman shock term and the full dispersion landscape 
Dint) and we keep the previously determined value of the 
Raman shock time (20 fs). The simulations are in good 
agreement and in particular, reproduce the experimen¬ 
tally observed cancellation and overcompensation of the 
soliton spectral recoil by the Raman induced red-shift. 

Conclusions. — In summary, we experimentally ob¬ 
served and studied for the first time the Raman-induced 
soliton self-frequency shift in SiN microresonators. Con¬ 
cerning the implications of our work, the results enable 
to increase substantially accuracy and understanding in 


numerical modeling of soliton based frequency combs in 
the soliton regime for SiN in particular, and amorphous 
microresonators in general. The observations should ap¬ 
ply to other classes of amorphous microresonators, such 
as aluminum nitride or silica based. Likewise the ob¬ 
servation of Raman induced red-shifted soliton spectra 
has implications for understanding the fundamental noise 
properties such as to the timing jitter of the emitted soli- 
tons (and thus phase noise in the generated microwave 
beatnotes), as the red-shift implies an energy loss to the 
phonon and results in additional noise as governed by the 
fluctuation dissipation theorem [44] . 
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